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Abstract We study the reactions 77 tt^tt", tt+tt", 
K^K^, K^K~ , rj and iv^rj based on a chiral Lagrangian 
with dynamical light vector mesons as formulated within 
the hadrogenesis conjecture. At present our chiral La- 
grangian contains 5 unknown parameters that are rel- 
evant for the photon fusion reactions. They parame- 
terize the strength of interaction terms involving two 
vector meson fields. These parameters are fitted to pho- 
ton fusion data 77 — >■ 1:^1:^, tt+tt", 77*^77 and to the de- 
cay r; — )• 7r'^77. In order to derive gauge invariant re- 
action amplitudes in the resonance region constraints 
from micro-causality and exact coupled-channel unitar- 
ity are used. Our results are in good agreement with the 
existing experimental data from threshold up to about 
0.9 GeV for the two-pion final states. The ao meson in 
the 7r°?7 channel is dynamically generated and an accu- 
rate reproduction of the 77 ir'^r] data is achieved up 
to 1.2 GeV. Based on our parameter sets we predict the 
77 — )■ K'^K'^, K~^K~ , rjr] cross sections. 

Keywords Meson production • Chiral Lagrangians • 
Vector-meson dominance • Partial-wave analysis • 
Dispersion relations 

PACS 13.60.Le • 12.39.Fe • 12.40.Vv • ll.80.Et • 
ll.55.Fv 

1 Introduction 

Photon- fusion reactions 77 — > PP (with PP = tt^tt'', 
TT+TT^, K^K'^, , rji] and n'^r]) play an important 

role in our understanding of non-perturbative QCD [1, 
. As a systematic approach chiral perturbation 
theory (xPT) is applied to describe these reactions at 
low energies [SJ[7] . Such studies have been performed at 
next-to-leading order (one loop) [SUH] and at next-to- 
next-to- leading order (two loops) (THHj- However, xPT 



is limited to the near-threshold region and cannot serve 
as an appropriate framework in the resonance region, 
where exact coupled-channel unitarity becomes an im- 
portant issue. In [TUllTTlfT^fT^ it was shown that the 
combination of electromagnetic gauge invariance, micro- 
causality and coupled-channel unitarity together with 
chiral symmetry helps to achieve a systematic approxi- 
mation of hadron interactions beyond the threshold re- 
gion. In the following the same approach is applied to 
the photon-fusion reactions. 

The cross sections of fusion processes are very sen- 
sitive to hadronic final-state interactions. Therefore, a 
crucial input for the present work is a proper descrip- 
tion of the reactions of two pseudoscalars into the same 
or different two pseudoscalars. Recently, within the novel 
unitarization approach developed in [lOlllllfTB] a coupled- 
channel computation of Goldstone-boson scattering has 
been performed ^2]. The calculations are based on the 
chiral Lagrangian formulated with light vector mesons. 

The light vector mesons play a crucial role in the 
hadrogenesis conjecture [T4 j [T5 l fT6 l fT7 l lT8 l [T9 l [20] . Together 
with the Goldstone bosons they are identified to be 
the "quasi-fundamental" hadronic degrees of freedom 
that are expected to generate the meson spectrum. For 
instance it was shown that the leading chiral interac- 
tion of Goldstone bosons with the light vector mesons 
generates an axial-vector meson spectrum that is quite 
close to the empirical one J,5j . It seems quite natural to 
keep vector mesons as explicit degrees of freedom when 
aiming for a description of hadron physics in the reso- 
nance region. In particular, this extends the resonance- 
saturation mechanism |21j by taking into account ex- 
plicitly the dynamics of the vector-meson propagator 
In the present analysis we further emphasize the 
significant role of light vector mesons in hadron physics. 
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One important aspect in photon-fusion reactions is 
the formation of scalar and tensor resonances. In partic- 
ular one finds the rather narrow scalar states ao(980) 
and /o(980) PSIITIIS]. In the present work we concen- 
trate on energies below 1.2 GeV. There the scalar res- 
onances /o(980) and ao(980) appear to dominate. In 
the approach of [H] both states are dynamically gener- 
ated from coupled-channel interactions, in accordance 
with the hadrogenesis conjecture. The two-pion chan- 
nel, where the /o(980) resonance is seen prominently, 
has been discussed in detail in IVF.. In contrast, the ttt] 
channel with the ao(980) state has not been addressed 
in [12]. There are no elastic ttt] scattering data avail- 
able. However, the tt rj channel can be populated by the 
inelastic photon-fusion reaction. We will show in the fol- 
lowing that the dynamical generation of the ao(980) is 
in very good agreement with the available experimental 
data on the photon-fusion production of tt rj. 

At higher energies tensor resonances come into play. 
In principle, also the tensor resonances /2(1270) and 
02(1320) are expected to be naturally generated within 
our approach from vector- vector interactions. However, 
wc focus in the present paper on lower energies and on 
the mutual interactions between Goldstone bosons. In 
the present approach vector mesons appear as exchange 
particles, but not as states in the coupled channels, 
i.e. we consider scattering and rescattering of the type 
77 — >■ PP and PP — PP, respectively, but disregard 
PP O VV, VP, where V generically denotes a vector- 
meson state. In principle, the basis for the systematic 
inclusion of vector mesons as coupled-channel states has 
been laid out in This technically more challenging 
computation is beyond the scope of the present work. 
Consequently the analysis of the resonances /2(1270) 
and 02(1320) and the corresponding energy regime is 
postponed to the future. 

Experimental information on the recations 77 — >■ 
hadrons is accessible in e+e" collisions via the reaction 
e+ + e- e+ + e'-h hadrons [25 l [26 l [27 l [28 l [29 l l30 l l3T] . 
Recently high-statistics data have been reported by the 
Belle Collaboration [5^[551[M1[55] . including the first 
measurement of 7777 production. The reaction 77 — >■ tt'' 77 
is linked to the decay rj — > 7r'^77 by crossing symmetry. 
Recent analyses of this decay have been performed at 
AGS [SSEZj and at MAMI [MIEH]- We will use the 
corresponding integrated and differential information 
about this decay to constrain our so-far unknown cou- 
pling constants parameterizing the strength of interac- 
tion terms with two vector meson fields. 

The main purpose of this paper is to obtain a unified 
description of the reactions 77 — )■ tt^tt", tt+tt", K^K^ , 
K'^ K~ , fjij and Tr'^ry using the dispersive effective field 
theory framework of [lOlfTTlfT^IT^ . The present study 



extends the analysis of Goldstone-boson scattering per- 
formed in [12]. In the next section we will specify our 
Lagrangian and comment on the determination of the 
corresponding coupling constants. Section [3] presents 
the required decomposition into partial- wave amplitudes. 
The unitary and causal summation scheme, which is at 
the heart of the coupled-channel dynamics, is reviewed 
in section [4] Numerical results and comparisons to data 
are given in section [5] and a summary of our work is pro- 
vided in section [6] 

2 Chiral interaction Lagrangian including 
vector mesons 

The study of the photon-fusion processes is performed 
in application of the chiral Lagrangian with explicit vec- 
tor mesons. The leading order terms of this Lagrangian 
have been constructed based on a formal power count- 
ing scheme |201I19| . The counting rests on the dynam- 
ical assumption of hadrogenesis and on large- A'^c argu- 
ments, where denotes the number of colors. In previ- 
ous successful applications two- and three-body decays 
of vector mesons have been presented in [T^IS^HO] . For 
photon-fusion reactions the relevant part of the leading- 
order Lagrangian takes the simple forirj^ 

C = -^^A^A^iv{<PQ[$, Q\_] 
+ i^^''tr{9^<Z' [Q, 
-efv^^A,tT{^^'''Q] 

+ '^d,A,tT{<P^'^[<P, [<P, Q]_]_} 

+ 'Jx^A^tr{[d,<P,<P^'^]_[Q,<P]_} (1) 

-Y^tr{5"<^."[['^^M-.<^''1J 

-3l72t'-{[^M.,M+[5"'^><^"1+} 

^Note that in | 12II19 | slightly different notations were used. 
The relations between ey, gn, 9f^ hp in |12II19| — denoted 
by old — and fv , gi, 92 , hp used here and in |20J are fv = 

0.776 GeV „ „ _ „ „ — „ or,^ !, — 0.776 GeV i, r„7 jl 
4^ ev, ffi — 9d, 92 — gp and hp — j- hp[old\. 
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where the Goldstonc-boson field <P, the vector-meson 

field ^/ij/, the charge matrix Q, and the mass matrix xo 

are normalized as follows: 
1 
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V2K+ 
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(2) 



Note that we assume perfect isospin symmetry through- 
out this work. In ([T]) the photon field is denoted by 
and e — 0.303 is the electromagnetic charge. For vec- 
tor mesons we use the antisymmetric tensor-field rep- 



resentation 



-^vii giving rise to the resonance- 



saturation mechanism |21| . In principle, in the hadro- 
genesis conjecture also the singlet eta field is part of 
the "quasi-fundamental" hadronic degrees of freedom. 
It can be included in the flavor matrix of in a straight- 
forward way [20 . However, as discussed in the introduc- 
tion, we do not include vector channels in our coupled- 
channel approach. Consequently we also do not include 
channels with the ry' which would appear in the same 
energy regime. Vector mesons are important, nonethe- 
less, as they contribute as exchange particles to the cou- 
pled channels 77 and PP. On the other hand, this is 
not the case for the rf . Consequently, here we have not 
considered the eta singlet explicitly in our Lagrangian. 

Finally we shall discuss the coupling constants ap- 
pearing in ([1]) . The following set of parameters has been 
determined by the masses and decay properties of the 
vector mesons [19,39 : 

/v = 0.140 ±0.014 GeV, /ia-2.10, 

/ip/y = 0.23GeV, 615 =0.92, (3) 

/~0.90GeV, 6a = 0.27. 

The values of the other parameters 51-3, 35 and ho 
have not been determined so far. Assuming that they 
are of natural size we will study in section [5] the impact 
of variations of these parameters on the photon-fusion 
processes and on the related decay 77 — >■ 7r"77. 



Table 1 The coupled-channel states characterized by 
isospin I and G-parity G. The Pauli matrices o-^ act on isospin- 
doublet fields K,K with for instance K = {K+,K°Y. Note 
that in particular the neutral (/s = 0) two-pion state with 
isospin two is given by -h= (2 Tr^ ttJ^ — -k^ -Kq — np tt^). 
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3 Partial-wave amplitudes 

We define the transition amplitude for the process 77 
PP as 



(4) 



(P(p)P(g)|T|A(fci,Ai)A(fc2,A2)) = 

{2-K)H\kr + k2-p~q) T^^" e^(fci, Ai) e.(fe, A2) 



where fci 2 and ei 2 are the momenta and the polariza- 
tion vectors of the incoming photons, respectively, and 
p, q are the momenta of the outgoing mesons. 

In general, the two-body scattering problem decou- 
ples into orthogonal channels specified by isospin, G- 
parity, parity and strangeness quantum numbers. For 
the case at hand, with two photons in the initial state 
and two pseudoscalars in the final state, parity is always 
positive (see below) and strangeness is always zero. In 
each of the channels, finally specified by isospin / and 
G-parity G, there are several meson-meson states cou- 
pled to each other. In Table [l] we have specified the 
states which contain the most relevant meson-meson in- 
formation below 1.2 GeV. Here we are neglecting multi- 
pion states which are only relevant for higher energies. 
In particular we neglect in this way also PV and VV 
states which on account of the resonant nature of the 
vector, V ^ states would significantly contribute to the 
multi-pion states. 

Lorentz covariance and gauge invariance lead to a 
decomposition of the scattering amplitude T^"^ into 
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Fig. 1 Tree-level diagrams for 77 — s> PP reactions with the exchange of pseudoscalar (solid line) and light vector (dashed line) 
mesons. 



Lorentz tensors L" and invariant amplitudes Fi, 

L't" = klk'i-{k,-k2)9>'\ (5) 
Lf = (Z\2 (fci ■k2)~2 {ki ■ A) {k2 ■ A)) .g^" 
-Z\2 k'( fc^ - 2(/ci • fcz) Zi^ A" 
+2(fc2 • /i) k'( A'' + 2(fci • A) fc^ Zi'' , 

where A — p — q and T'^'^ satisfies the Ward identities 



(6) 



The motivation for choosing these particular Lorentz 
structures is twofold. First, the corresponding invariant 
amplitudes are independent and free of kinematical sin- 
gularities or zeros. Second, in order to simplify further 
calculations we have chosen the Lorentz tensors such 
that the following property holds: 



2 fiu 



(7) 



The invariant amplitudes Fi and F2 are analytic 
functions of s, t and u except for dynamical cuts. Fur- 
thermore, a t — li crossing symmetry is satisfied due to 
the Bose statistics of the two photons. 

It is useful to introduce the helicity components of 
the scattering amplitude and decompose each of them 
into their partial waves, 

<^++ =r^''e^(fci,+l)e,(fc2,+l) 

= J2 (2J+l)ii'i4o'^(cos0), 

even J>0 

-r''''e^(fci,+l)e,(fc2,-l) 

= J2 i2J+l)t^^ldg\cos0), (8) 

even J>2 



where 



e^(ki,A = ±l) 



/ \ 



_ I 

V J 



e,(k2 = -ki,A = ±l) = 



and d{^j^{cos 9) are Wigner rotation functions. In (|8l) 9 is 
the center-of-mass scattering angle. The axis of tne col- 
liding photons in their center-of-mass frame has been 
chosen to agree with the third axis of the coordinate 
system. Note that the partial-wave expansion involves 
only even J > A and positive parity P = -\- [41]. This 
constraint arises from the combination of Bose symme- 
try of the two initial massless photons with the possible 
J, P quantum numbers of the two Goldstone bosons in 
the final state. 

The invariant amplitudes Fi^2 can be expressed in 
terms of the helicity amplitudes 4'++i 4>-\ i 



/ \ 

I 

V J 









1 



(9) 



where Pcm is the final center-of-mass relative momen- 
tum and X = cos 9. For unpolarized photons the differ- 
ential cross section is given by 



dcos9 32 7rs 4 



(10) 



where /3 — 2pcm/\/s- If two identical particles appear 
in the final state (neutral pions, for instance) one has 



to include an additional factor of 1/2 in ( 10 ) or perform 
the integration only over 6* S [0, 7r/2]. 

According to ([8| the partial-wave helicity ampli- 
tudes t^-'j, 4- 

can be computed in terms of the in- 
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variant amplitudes Fi_2 as 



dx 



sFi(s,a;)4o^(x), 



dx 



2pl„,s {x^-l) F2{s,x)dg\x), (11) 



with the help of the useful identities for the Wigner 
rotation functions HH 



d\^^{x) = Pj{x)., 

dii\x) = 



2xP'j{x) 



v/(J-l) J(J+1) (J + 2) 



(J-l)(J + 2) 



(12) 



In order to avoid kinematical singularities and zeros 
in the partial-wave amplitudes at threshold, we rescale 
(11) by a phase-space factor [pcmPcmY , 



T 



(./) 



iJ) 



(.Pcni Pcm , 



J ++ ' 



yycrn ycm J 

and also multiply by s"^ to ensure a finite limit of the 
phase-space matrices at large energy. 

The invariant amplitudes (|5| computed from the 
chiral Lagrangian ([T]) read 
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2 rri! 
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(14) 



where the sum runs over the octet of Goldstone bosons 
([8]) or the vector-meson nonet ([9]) and m^.y^z de- 
notes their respective masses. The coefficients C... are 
presented in Tables [3) |4j [5) and [6] with respect to the 
coupled-channel states of Table [T] In Fig. [T] the set of 
tree-level diagrams that gives nonzero contributions is 
depicted. Note that for the isospin states which contain 
identical particles (e.g. Itttt), I?/??), . . . ) we use a con- 
vention where the unitarity condition for identical and 
non-identical two-particle states are the same. 

The conventions of Table [l] imply the following rela- 
tions between scattering amplitudes in isospin and par- 
ticle bases, 



T, 



-,1=0 



77— ^TT^TTT 



= 2 f — T' 



Tll = 2 

77— >7r TV I ' 



-^77— ^tt'^tt 



F 



77-i. K+K- 



T. 



1=0 

77-> K K ^ 77->- K K ' 



T. 



2 rj,I = 2 
2 77— VTT TT 



1 77^. jjr] — ^ -L 77-> r; ?; i 



rri ^ rpI — rpI—\ 



(15) 



where the factor 2 reflects our normalization for two- 
body states with identical particles. Note that in [T^] 
meson-meson interactions were studied only for S- and 
P-waves. The result for the D-wave amplitude can easily 
be obtained from Eq. (4) and Eq. (5) of [12 . 



The partial- wave amplitudes obtained from (13 11 ) 



and ( 14 1 at tree-level will serve as an input for the non- 
perturbative coupled-channel calculations to which we 
turn next. 
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Table 2 The positions of the closest left-hand branch points 
of t^outside(s) for 77 — > J'f that are determined by the t- and 
n-channel exchange processes. The numbers in the column 
"ch = ab" correspond to the out-states (a) and in-states (b) 
of Table [U 



Fig. 2 Locations of left-hand cuts of the 77 -> tttt partial- 
wave amplitude in the complex s-plane. The branch point of 
the p- meson exchange is located at = — (m^ — rn^)-^ /m"^, 



while one- and two-pion exchange cuts start at 



0, 



— 9m^/4, respectively. The dashed line identifies the conver- 
gence region of the conformal expansion \20\ . 

4 Coupled-channel dynamics 

The partial-wave reaction amplitudes satisfy the dispersion- 
integral representation, 

TL{s) = Ui{s) (16) 





ch. 






Description 
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ml 
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t, M-ch (27r) 
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31 
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t, M-ch (27r) 

t, M-ch {ttK) 



2+ 



21 



9 m-^ 

4r 



t, M-ch (27r) 



E 



where the phase-space matrix p'^^{s) is diagonal in c and 



d. In ( 16 ) the coupled-channel indices a and b run over 
the various channel 77, tttt, etc. The unitarity con- 
dition implies that the discontinuity of the generalized 
potential 



2i 

= for 



S > Mthr 



vanishes for energies larger than the s-channel thresh- 
old. This is nothing but the condition that C//f,(s) has 
left-hand cuts only. 

In our normalization the hadronic part of the phase- 
space matrix 



1 

8^ 



Pcm 



2 J+1 



(18) 



approaches a finite value in the high-energy limit. Inter- 
mediate states with two photons are neglected in this 
work. Numerically they are largely suppressed being 
proportional to at least. Therefore the two-dimensional 
phase-space matrix for the two-photon states need not 
to be specified here. 

Our approach satisfies the electromagnetic gauge- 
invariance constraint. This follows from the on-shell 
condition for the generalized potential, for which we will 
construct a systematic approximation in the following. 
The on-shell reaction amplitude will then be derived in 



application of ( 16 ). Owing to the matching scale pM in 



( 16 ) the non-perturbative coupled-channel calculation 



chiral Lagrangian smoothly connect at s = pfj. This 
is discussed in detail in [TU] and [T^. Here we iden- 
tify fiM with the smallest two-body hadronic threshold 
value and thereby assume the applicability of xPT at 

s = fllj. 

Following TIT the generalized potential U^ij{s) can 
be extrapolated to higher energies in a controlled man- 
ner by applying conformal mapping techniques. In a 
first step the generalized potential is split into two con- 
tributions 



(17) U{s) = C/insidc(s) + [/outsido(s) , 



(19) 



where J7insidc(s) contains the contributions from close- 
by left-hand cuts and C/outsidcls) the contributions from 
far-distant left-hand cuts. While the former can be ex- 
plicitly calculated from the chiral Lagrangian in a per- 
turbative application, the latter reflect short distance 
physics that need to be parameterized systematically 
and efficiently. For the specific example reaction 77 — ^ 
TT TT the separation of 'inside' and 'outside' is implied by 
the dashed line in Fig. [2j 

The outside potential is expanded in powers of a 
conformal variable ^(s) constructed such as to ensure 
convergence for any value of s inside the area bounded 
by the dashed line of Fig. [2| 



Uo 



itsidc 



is) 



k=0 



Ck^Hs) for s<Al, 



(20) 



where the coefhcients Ck are uniquely determined by the 
first k derivatives of L'outsidc(s) at the expansion point 
s = fi^ . We identify this expansion point with the mean 



^In the following we do not display any more the angular- 
momentum superscript J explicitly. Where not needed we 



and the results from a perturbative application of the also do not display the channel index ab. 
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of the initial and the final threshold. Following pQ] the 
function ^(s) takes the form 

a{Al-sr-l 



(a-2 6)(yl2-s)2 + l 



(21) 



where the parameter Aq is determined by the posi- 
tions of the closest left-hand branch point of the struc- 
tures which enter J/outsidc- The quantities iie (expan- 
sion point for outside potential) are collected in Table 
[2] for the isospin states of Table [T] Finally, As is the 
upper limit of the convergence region, see Fig. [2] 

To be specific, the 'inside' part of the potential re- 
ceives contributions from the one-pion (kaon) exchange 
processes only. We evaluate the contributions from the 
cuts starting from 

9to2 



< s < 



for 



77 



2 771^)2 



< s < for 77 -> KK , (22) 



2 {"^TT 

— TTl 

(m^ + niK) 

following the procedure outlined in Appendix B of |10j . 
According to Table [3] there are no more cases to be 
considered. The coefficients Ck in the outside part of the 
potential are computed by evaluating the first n deriva- 
tives of the partial-wave amplitudes as determined via 



(11 12) by the tree-level resuh (14) 



Via ( 19 ) we obtain an approximated generalized po- 
tential for energies Aq < s < A^. While the inside part 
of the potential is defined for s > A^ also, the outside 
part is undefined for s > A^ by (20). We may simply 
cut off the integral in ( 16 ) at s = However, it is ad- 
vantageous not to do so, since that would induce rapid 
variations of the amplitudes close to and below s = Ag. 
While the precise form of the generalized potential at 
s > should not influence the reaction amplitudes 
in the target region, where we have a controlled expan- 
sion, it is useful to minimize its residual influence on the 
target region. This is the case if the outside potential 
is smoothly extended for s > vl^ by a constant (see [TUl 
[TTIIT^IT^ ). We note that due to the particular form of 



the conformal map (21), the generalized potential and 
its derivative are continuous at s = A^. 

Following [121113]) we choose in the following 71 = 3 
in (20 1 and As = 1.6 GeV. We have varied As in the 



range 1.4 GeV to 1.8 GeV. The impact of this variation 
on the results is small [43] . 

A crucial observation behind our summation scheme 
is the fact that the computation of the partial-wave 



scattering amplitude from ( 16 ) at energies larger than 



threshold also requires only the knowledge of the gener- 
alized potential at energies larger than threshold. More 
generally, depending on where we want to compute the 



Table 3 The coefficients Csg, C^""^ and C^^' = 2 C^^' of the 
invariant amplitudes ( |14| ) with respect to the coupled-channel 
states of Table [ij The numbers in the column "ch = ab" 
correspond to the out-states (a) and in-states (b) of Table [ij 
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partial-wave scattering amplitudes, it suffices to con- 
struct a controlled approximation of the generalized po- 
tential in a specific region of the complex plane only. 



This is always achieved with (19 20 ) and the desired 



solution of ( [16| can be found by the N/D ansatz [H] 
T,fc(s) = ^Z?-i(s)iV,b(s), (23) 

c 

where Dab{s) contains only the right-hand s-channcl 
unitarity cuts, 

ds s- Nac{s)pcb{s) 



Dab{s) 



E 



(24) 



and the matrix Nab{s) contains left-hand cuts only, 

Nab{s) = Uab{s) (25) 



E 



d± s- iilj Udbjs) - Udbjs) 



The system ( 25 ) with the input ( 19 ) , ( 20 1 can be solved 



numerically by the method of matrix inversion. 



5 Numerical results 

In this section we present our results for the cross sec- 
tionfj^of the reactions 77 7r°7r°, tt+tt^, K^K'^, K+K~ , 
•q rj and tt^t;, evaluated with the J = 0, 2 partial- waves 
amplitudes. We have checked that the contributions 
from the higher partial waves are negligible in the en- 
ergy range ^/s < 1.2 GeV. 

^Usually the experimental results are limited to a range of 
|x| < Z with X = cos 9. In this case the cross section is given 



by a = 2 / 




■ dx. 
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We use the set of parameters given in (|3| for all the 
numerical results. However, the remaining five param- 
eters gi, 52, 33, <?5 and ho have to be determined. Our 
strategy is to use the empirical data on the reactions 
77 — >■ 7r°7r°, tt+tt" and TT^rj and in addition the dif- 
ferential and integrated data for the decay -q — > 7r°77. 
The results for the cross sections of the reactions 77 — > 
K^K^, K^K~ and 7777 are then pure predictions. 

On account of crossing symmetry, the decay ampli- 
tude for 77 — >■ 7r'^77 can be easily obtained from 77 — > 
Tr'^ry by considering tt*^ and photons as outgoing parti- 
cles. To get the decay amplitude it is enough to replace 



(fci 

{P- 
(P- 



+ k^f (fci 



yk^Y 



77 ' 

-'"71 7r , 

72"' ' 



(26) 



in the invariant amplitudes ( 14 ) . The differential decay 
rate is given by [45j 



dr = 



1 



1 



(27r)3 32 m: 



,3 

pol 



(27) 



For the integrated partial decay width one has to in- 
clude the degeneracy factor of 1/2 to account for the 
fact that one has two indistinguishable photons in the fi- 
nal state. To obtain the decay amplitude we use directly 



the tree- level result ( 14 ) . Since we are here in the low- 



energy decay region, we assume that coupled-channel 
effects are less important. For the reaction amplitudes 
we use, of course, the full rescattering formalism out- 
lined in the previous section. 



In a first step we use the reaction data of 77 — 
vr^TT^, 7r"'"7r~ and TT^r] to correlate the five free parame- 
ters. Having matched the data with the coupled-channel 
calculations leads to the following relations: 



gi = 0.900 - 0.200 53 + 0.038 hf^ + 0.128 ho , 



32 



1.50-0.2733 + 0.2535. 



(28) 



This leaves us with three free parameters. If they are 
varied within the range 173, 55, ho S [—5, 5] one obtains 
the cross sections depicted in Fig. |3] A detailed discus- 
sion of the cross sections will be given below. To get a 
feeling for the influence of the five parameters we also 
provide the cross sections for the case where all these 
five parameters are put to zero, see the dashed lines 
in Fig. [3] Obviously, one would significantly underesti- 
mate the data in both neutral channels tt'^tt" and Tr'^jy 
without the parameters gi and ho- Note, however, that 
the qualitative structure does not depend so much on 
these parameters. 

We continue with a determination of the remaining 
parameters using the existing data on ry — )■ 7r'^77 decay. 
The present experimental status for r] — >■ 7r°77 decay is 
the following: The Particle Data Group [IS] gives the 
branching ratio / — (2.7 ± 0.5) • 10~* and 

the full width = (1.30 ± 0.07) keV. This results in 
a partial decay width of r'^_>.^o^^ « (0.35 ± 0.09) eV. 
Theoretical studies have been performed in [T7lH51H5] . 

For the decay rj — > 77^77 three of the yet undeter- 
mined parameters contribute, namely 31, 33 and ho- 



Using the relation (28) for gi, we adjust 33 and ho to 
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the partial decay width and to the two-photon invariant- 
mass distribution depicted in Fig.|4j In this way we find 
33 = -4.88 , ho = 3.27 (29) 
which imphes gi = 2.70 and ,92 = -0.18 + 0.25 55. The 
fit yields 

r^^^o^^ = 0.310 eV (30) 

for the integrated partial decay width, in good agree- 
ment with the experimental value. 

We have determined four of our five free parameters. 
In the following we will show results where the remain- 
ing free parameter 35 is varied in the range (75 G [—5, 5]. 



Note that the achieved determination of the parame- 
ters is also crucial for future investigations. Originally 
all these parameters concern interactions between two 
vector mesons and an odd {hp) or even {gi, i ^ 1, 2, 3, 5) 
number of Goldstone bosons. In the future it is planned 
to explore also the importance of vector-meson channels 
for the coupled-channel problems (cf. the corresponding 
discussion in the introduction and [23^ ). There the cou- 
pling constants hp and gi enter directly and mediate, 
e.g., the transition from two vector to two pseudoscalar 
mesons. For the following reason these coupling con- 
stants are also important for our case at hand, in spite 



10 



Table 6 


The coefficients C^f'^\ C^^'^^ 


> ^93 > ^95 


and C<"'^'. 
On 


See the caption of Table | 


^for more 


details 






ch. 


(P.P) 




(P,0) 


(a;, uj) 


{4>,uj) 




(0,</') 




0+ 


21 
31 
41 


16 

8 

16/3 














16 

8 
9 

16 

27 




16 

"9" 








32 
9 
128 
27 


1- 


21 





16 [2 
3 V 3 


















31 





8 
3 


16 

3 













2+ 


21 





















'-'a-i 


0+ 


21 


32 





















31 


8 








8 

9 


16 
9 




32 
9 


1" 


21 





8 
3 


16 

3 













2+ 


21 























0+ 


21 


4v^ 








4 
3^/3 












31 


4 








4 
9 







16 
9 




41 


4 
3 








4 

27 







32 
27 


1- 


21 





4 [2 
3 V 3 


















31 





4 
3 
















'-■95 


0+ 


21 


8 

-2 





















31 








2 
9 


4 

9 




8 
9 


1- 


31 





2 
3 


4 
3 













2+ 


21 























0+ 


21 


32x/3m2 








3^/3 -^ 












31 
41 


32 m|, 

3 T 










9 

27 T 






32 
27 


9 ™if 
(16m2, -8m2) 


1- 


21 





32 /i" 2 


D 


D 












31 





3 '"if 

















11 



300 



I 


Markll 


I 


Belle 


I 


CELLO 




Belle syst. eiTor 



yy 



71 71 




60 



40- 



20- 





Crystal Ball 




BeUe 




Belle syst. error 




Fig. 3 Total cross sections for 77 — s> tt+tt witli 
|cos6»| < 0.6 (top left), 77 -s- tt°it° with |cos6»| < 0.8 
(top right) and 77 — ^ tt"?? with |cosS| < 0.9 (bot- 
tom). A variation of parameters gs, gz, ho S [— 5i 5] 
using \2fi\ is reflected by the various bands. Setting 
Qi = = ho yields the dashed lines. The data are 
taken from | 57ll3^ll25ll56ll33liaSII54 |. 
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Fig. 4 The single-differential invariant-mass distribution of 
the decay r) —^ ^"77. Parameters are chosen according to ( |28| , 
( |29[ |. Note that the parameters g2 and g^ do not contribute 
to this decay. The data are taken from |36ll37ll46l . 



of the fact that we do not consider the vector-meson 
channels: The neutral vector mesons couple directly to 
photons; see also Fig. [l] Therefore, the coupling con- 



stants ho and gi enter also the transition amplitudes 
from two photons to two pseudoscalars. In turn, data 
on such interactions between hadrons and electromag- 
netism can be used to constrain purely hadronic cou- 
pling constants. This resembles our determination of Ha 
from the decay to — 771" in ^9j. Note, however, that in 
our formalism this line of reasoning does not lead auto- 
matically to strict vector-meson dominance, but rather 
to an improved version thereof [iUll^U] . 

We now turn to a detailed discussion of the vari- 
ous two-meson channels populated by photon fusion. 
The first highlight is the ir'^rj channel depicted in Fig. 
[5| Here our formalism shows a dynamically generated 
scalar-isovector resonance which is in full quantitative 
agreement with the experimental data; see also [50) 
I23ll51j where similar findings have been reported. In 
our approach we find that this ao(980) resonance co- 
incides with the two-kaon threshold and emerges from 
rescattering and coupled-channel effects between 77^77 
and KK. We contrast our full coupled-channel result 



with a pure tree-level calculation based on (14 1. The 



12 




60 



40- 



20- 



+ 


Crystal Ball 


i 


BeUe 




Belle syst. error 



s'"[Ge\^ 




s'"[Ge\^ 



Fig. 5 Cross section for the reaction 77 — ^ Tr'^r; using ( |28[ ), ( |29[ l togetlier with 35 g [—5, 5]. The full result is shown on the left 
and the tree-level result on the right. The tree- level result does not depend on 35. Also the dependence of the full result on 
is very weak. It is caused by the cross-channel effect 77 — 5> KK — !> tt^'tj. See the figure caption of Fig. [3] for more details. 



latter is also depicted in Fig. [5] and, of course, does not 
show a resonance shape, in obvious disagreement with 
the experimental data. We stress again that according 
to the hadrogenesis conjecture [14, 15, 16, 17, 18lll9l[20| 
the low-lying scalar resonances are supposed to be gen- 
erated dynamically. An incarnation of this proposition 
is seen in Fig. [5] We recall from our previous discussion 
about Fig.[3]that the location of the resonance does not 
depend on the choice of the coupling constants gi and 
ho- Only the height of the curve is sensitive to these 
parameters. This provides confidence in the robustness 
of our interpretation of the lowest-lying scalar-isovector 
resonance. 

The cross sections for the two-pion channels are de- 
picted in Fig. [6] Obviously both channels tt+tt" and 
tt^tt" are well described up to energies of about ^/s w 
0.9 GeV. Then our calculations show a distinct peak, 
most pronounced in the neutral channel. After this peak 
our theory curves decrease while the data continue to 
rise. Two issues need to be disentangled here, namely 
the location of the /o(980) in the S-wave and the rise 
towards the tensor mesons in the D-wave. To do this we 
compare our results also to the partial-wave analysis of 
[I] as shown in Fig. [t] For the D-waves (bottom panels) 
we observe reasonable agreement up the point in energy 
where the peak from the isoscalar tensor meson starts 
out. As already stressed in the introduction we expect 
that in the spirit of the hadrogenesis conjecture this 
peak will be generated by vector-vector channels. But 
since this is beyond the present work we cannot expect 
to obtain a reasonable description of the D-wave beyond 
about 0.9 GeV. Below this energy the agreement is very 
satisfying. Turning to the S-wave we observe also good 
agreement for isospin 1 — 2 (top right panel in Fig. 



[7|. For the isoscalar channel (top left) some disagree- 
ment with the results of [1] is observed. Most notably 
our peak for the /o(980) is slightly shifted to lower en- 
ergies, i.e. this dynamically generated scalar-isoscalar 
state is somewhat overbound in our approach. This has 
already been observed in il2 ■ Whether this is due to 
higher-order effects in the scattering kernel or due to 
missing vector-vector channels remains to be seen. 

For the energy range below 0.9 GeV we deduce from 
Fig.[6]that we have obtained an overall good description 
of the reaction data. This is completely in line with the 
complementary information contained in the pion phase 
shifts as addressed in [T^]. In Fig. [s] we compare our cal- 
culations to the results from chiral perturbation theory 
(xPT)[8l|9l[7l|6|. We observe satisfying agreement. Note 
that even without vector mesons our calculations con- 
tain multi-loop diagrams by the achieved resummation 
in the s-channel. On the other hand, our calculation 
does not contain all one-loop diagrams in the t- and 
u-channel which enter at next-to-leading order. 

In view of these differences one can be satisfied with 
the agreement and conclude that the numerically most 
important corrections from the xPT point of view are 
included in our approach. At larger energies pure xPT 
ceases to work and resummations must be incorporated 
in one or the other way [5U 1 E5 1 [5 ^ [T ^ H5] . 

Finally we show in Fig. |9] the result of a tree-level 
calculation based on our amplitudes ( 14 1. Obviously the 



charged-pion channel is fairly insensitive to rescattering 
effects, i.e. to a large extent dominated by the one-pion 
exchange, which is responsible for the steep rise of the 
cross section at low energies. The neutral-pion chan- 
nel, however, which does not have the corresponding 
one-pion exchange, is dominated by loop/rescattering 
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Fig. 6 Cross section for the reactions 77 — > tt+tt" (left) and 77 — > tt^tt" (right) using ( |28[ ), ( |29[ l together with 35 G [—5,5]. 
The obtained region is limited by the full thick lines. For the charged pions (left) the S-wave (long-dashed) and D-wave 
(short-dashed) are shown separately. See the figure caption of Fig. |3] for more details. 
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Fig. 8 Comparison of our results (full lines) to the calculations from chiral perturbation theory (xPT) and to data. The dashed 
lines denote the next- to-next-to- leading-order calculations. The dotted lines denote the respective lowest-order non-trivial xPT 
result, which is leading order for the charged case and next-to-leading order for the neutral case. See also the figure caption of 
Fig. [sjfor more details. 
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Fig. 9 Comparison of tree-level calculations to the full results and to data. Tree-level calculations are depicted by dashed 
lines. See also the figure caption of Fig. |3] for more details. 



effects. In this channel the tree-level calculation fails 
already at low energies. 

While the previously discussed channels have been 
used to some extent to fix our free parameters, the chan- 
nels which we discuss in the following are pure predic- 
tions. Unfortunately the data situation is rather poor 
in all three channels 77 — K~^K~, K^K'^ and rjrj, but 
we will see that it is a non-trivial task to match the 
available data points. We restrict ourselves to the en- 
ergy region close to threshold, i.e. to ^/s < 1.2 GeV. 
There we expect the S-wave to dominate such that we 
do need to worry about the tensor mesons. 

The reaction 77 — >■ K~^K~ is depicted in the top 
left panel of Fig.jlOj Other theory approaches have been 
reported in [23, 53, 54j . Unfortunately there is only one 



data point with a large energy uncertainty in the con- 
sidered energy interval. Nonetheless, this data point is 
significantly lower than generic tree-level calculations. 
For comparison we show two types of such tree-level cal- 
culations. The dashed line is obtained if our Lagrangian 
is used directly for the amplitude and not for the poten- 
tial of the full coupled-channel calculation. An alterna- 
tive tree-level approach is to use just the kaon-exchange 
Born diagrams. We recall that the corresponding pion- 
exchange Born diagrams are very significant for the low- 
energy part of the reaction 77 tt+tt^. (This is the 
XPT-LO curve of Fig.jsj) For the kaon case the situation 
is obviously different. While tree- level calculations fail 
to reproduce even the close-to-threshold data, our full 
coupled-channel approach leads to a significant reduc- 
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Fig. 10 Cross sections for the reactions 77 — s> K'^ K~ 
(top left), 77 — > (t op right) and 77 — > rjjj (bot- 
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The obtained region is limited by the full lines. Tree- 
level calculations are depicted by dashed lines. For the 
charged-kaon case a pure Born-term calculation (one- 
kaon exchange) is shown by the dotted curve. The data 
are taken from I29II30||5T]|35I . 
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tion of the Born amplitude and matches the available 
data point very nicely. Hence, the final-state interac- 
tions are strong in this channel. A similar finding has 
been reported in [23 . Finally we note that our approach 
shows a visible cusp at the two-eta threshold. It is even 
more pronounced in the neutral-kaon channel to which 
we turn next. 

The reaction 77 ^ K^K^ is shown in the top right 
panel of Fig. [lO] The data, albeit with large error bars, 
point to an initial steep increase of the cross section 
with energy followed by a not so rapid fall. This behav- 
ior is qualitatively reproduced by our full calculation, 
though we do not fully match the second data point 
quantitatively. Tree-level calculations cannot reproduce 
at all this rise-and-fall behavior. Indeed, it is natural to 
expect that final-state interactions are strong because 
the photons couple stronger to the intermediate charged 
states than to the final neutral ones. 

The bottom panel of Fig. [TO] shows the cross sec- 
tion for the reaction 77 — )• ryry (see also [SI]). The data 
suggest a rather flat energy dependence which cannot 



be reproduced by a pure tree-level calculation (dashed 
line). In contrast, our full calculation (full lines) includ- 
ing rcscattering meets this requirement of a compara- 
tively flat cross section. 

We stress again that our results for the two-eta and 
two-kaon channels are pure predictions. These chan- 
nels did not enter the determination of free parame- 
ters. Note also that the results are basically insensitive 
to the remaining free parameter g^. Clearly, better data 
in these channels would be highly welcome to further 
check the validity of our coupled-channel approach with 
dynamical vector mesons. 



6 Summary and outlook 

We have performed a controlled study of the reactions 
77 — 7r°7r°, 7r+7r~, K^K^, K^K^, rj-q and ir^rj in the 
energy regime between the respective thresholds and 
about 1.2 GeV. The reaction amplitudes were derived 
from the chiral Lagrangian with dynamical vector me- 
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son fields properly constrained by micro-causality and 
coupled-channel unitarity. 

There are 5 unknown parameters, which have been 
constrained from the reactions 77 — >■ tt'^tt'^, Tr~^iT~ , ir'^r] 
and from the differential decay t] — >■ 7r°77. In particu- 
lar we have achieved an excellent description of the re- 
action 77 — )• TT^r] with its lowest-lying scalar-isovector 
ao(980) resonance. The ao(980) resonance peak position 
does not depend on any of the free parameters. Based 
on our parameter constraints we predict the low-energy 
77 — > K^K~ , K^K^ and 7777 cross sections. 

While the vector mesons do play a crucial role in the 
derivation of the generalized potentials for 77 — )• PP 
and PP ^ PP with P = TT,K.ru the feedback of 
PP -J> PV,VV reactions with V = p,uj,K*,(j> re- 
mains to be studied systematically. According to the 
hadrogcncsis conjecture we expect a quantitative de- 
scription of the reactions 77 7r°7r°, tt+tt", K^K^, 
K^K~ , T] 7] and tt^t] up to about 2 GeV once such chan- 
nels are incorporated in a controlled manner. 

Accurate low-energy photon fusion data in particu- 
lar in the strangeness channels would further scrutinize 
the intricate three-flavour dynamics of the Goldstone 
bosons and light vector mesons. 

Acknowledgements The work of CT and SL has been sup- 
ported by the European Community Research Infrastructure 
Integrating Activity "Study of Strongly Interacting Matter" 
(HadronPhysicsS, Grant Agreement No. 283286) under the 
Seventh Framework Programme of the EU. 



References 

1. M.R. Pennington, T. Mori, S. Uehara, Y. Watanabe, Eur. 
Phys. J. C 56, 1 

2. J. A. Oiler, L. Roca, Eur. Phys. J. A 37, 15 

3. R. Garcia-Martin, B. Moussallam, Eur. Phys. J. C 70, 
155 

4. M. Hoferichter, D.R. Phillips, C. Schat, Eur. Phys. ,J. C 
71, 1743 

5. Y. Mao, X.G. Wang, O. Zhang, H. Zheng, Z.Y. Zhou, 
Phys. Rev. D 79, 116008 

6. J. Gasser, M.A. Ivanov, M.E. Sainio, Nucl. Phys. B 745, 
84 (2006) 

7. J. Gasser, M.A. Ivanov, M.E. Sainio, Nucl. Phys. B 728, 
31 

8. J. Bijnens, F. Cornet, Nucl. Phys. B 296, 557 

9. J.F. Donoghue, B.R. Holstein, Y.C. Lin, Phys. Rev. D 
37, 2423 

10. A. Gasparyan, M.F.M. Lutz, Nucl. Phys. A 848, 126 

11. I.V. Danilkin, A.M. Gasparyan, M.F.M. Lutz, Phys. Lett. 
B 697, 147 

12. I.V. Danilkin, L.I.R. Gil, M.F.M. Lutz, Phys. Lett. B 
703, 504 

13. A.M. Gasparyan, M.F.M. Lutz, B. Pasquini, Nucl. Phys. 
A 866, 79 

14. M.F.M. Lutz, E.E. Kolomeitsev, Found. Phys. 31, 1671 
(2001) 



15. M.F.M. Lutz, E.E. Kolomeitsev, Nucl. Phys. A 730, 392 

16. M.F.M. Lutz, E.E. Kolomeitsev, C.L. Korpa, Prog. 
Theor. Phys. Suppl. 156, 51 (2004) 

17. M.F.M. Lutz, E.E. Kolomeitsev, Nucl. Phys. A755, 29 
(2005) 

18. M.F.M. Lutz, M. Soyeur, Nucl. Phys. A813, 14 (2008) 

19. M.F.M. Lutz, S. Leupold, Nucl. Phys. A 813, 96 

20. C. Terschliisen, S. Leupold, M.F.M. Lutz, arXiv:1204. 
4125 [hep-ph] 

21. G. Ecker, J. Gasser, A. Pich, E. de Rafael, Nucl. Phys. 
B 321, 311 

22. B. Kubis, U.G. Meissner, Nucl. Phys. A 679, 698 

23. J. A. Oiler, E. Oset, Nucl. Phys. A 629, 739 

24. M.F.M. Lutz, I. Vidana, Eur. Phys. J. A 48, 124 

25. J. Boyer, et al., Phys. Rev. D 42, 1350 

26. H. Behrend, et al., CELLO Collaboration, Z. Phys. C 56, 
381 

27. H. Marsiske, et al.. Crystal Ball Collaboration, Phys. 
Rev. D 41, 3324 

28. D. Antreasyan, et al.. Crystal Ball Collaboration, Phys. 
Rev. D 33, 1847 

29. H. Albrecht, et al., ARGUS Collaboration, Z. Phys. C 
48, 183 

30. H.J. Behrend, et al., CELLO Collaboration, Z. Phys. C 
43, 91 

31. M. Althoff, et al., TASSO Collaboration, Z. Phys. C 29, 
189 

32. S. Uehara, et al.. Belle Collaboration, Phys. Rev. D 79, 

052009 

33. T. Mori, et al.. Belle Collaboration, J. Phys. Soc. Jap. 

76, 074102 

34. S. Uehara, et al.. Belle Collaboration, Phys. Rev. D 80, 
032001 

35. S. Uehara, et al.. Belle Collaboration, Phys. Rev. D 82, 
114031 

36. S. Prakhov, eConf C070910, 159 (2007) 

37. S. Prakhov, et al., Phys. Rev. C 78, 015206 

38. M. Unverzagt, Crystal Ball at MAMI Collaboration, 
Nucl. Phys. Proc. Suppl. 198, 174 

39. S. Leupold, M.F.M. Lutz, Eur. Phys. J. A 39, 205 

40. C. Terschlusen, S. Leupold, Phys. Lett. B 691, 191 

41. V.B. Berestetskii, E. Lifshitz, L.P. Pitaevskii, Quantum 
Electrodynamics (Pergamon Press, 1982) 

42. D. Varshalovich, A. Moskaev, V. Khersonskii, Qunatum 
Theory of Angular Momentum (World Scientific, Singa- 
pore, 1988) 

43. I.V. Danilkin, M.F.M. Lutz, arXiv:1208.2568 [hep-ph] 

44. G.F. Chew, S. Mandelstam, Phys. Rev. 119, 467 

45. K. Nakamura, et al., J. Phys. G G 37, 075021 

46. K. Lalwaiii. Measurement of the branching ratio of a rare 
decay rj n""/"/ with WASA-at-COSY. Ph.D. thesis. 
Department of Physics, Indian Institute of Technology 
Bombay (2010) 

47. L. AmetUer, J. Bijnens, A. Bramon, F. Cornet, Phys. 
Lett. B 276, 185 

48. E. Oset, J.R. Pelaez, L. Roca, Phys. Rev. D 67, 073013 

49. E. Oset, J.R. Pelaez, L. Roca, Phys. Rev. D 77, 073001 

50. J. A. Oiler, E. Oset, Nucl. Phys. A 620, 438 

51. M. Doring, U.G. Meissner, E. Oset, A. Rusetsky, Eur. 
Phys. J. A 47, 139 

52. A. Gomez Nicola, J.R. Pelaez, Phys. Rev. D 65, 054009 

53. N.N. Achasov, G.N. Shestakov, Phys. Usp. 54, 799 

54. C.H. Lee, H. Yamagishi, I. Zahed, Nucl. Phys. A 653, 
185 



